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ABSTRACT 

The determination of the Ba abundance in globular cluster (GC) stars is a very powerful test to 
address several issues in the framework of multiple population scenarios. We measured the Ba content 
for a sample of more than 1200 stars in 15 Galactic GCs, using high-resolution FLAMES/Giraffe spec- 
tra. We found no variation in [Ba/Fe] ratios for different stellar populations within each cluster; this 
means that low-mass asymptotic giant branch stars do not significantly contribute to the intracluster 
pollution. 

Very interestingly, we obtained that the fraction of Ba stars in first generation (FG) stars is close 
to the values derived for field stars (~2%); on the other hand, second generation (SG) stars present 
a significant lower fraction. An independent and successful test, based on radial velocity variations 
among giant stars in NGC 6121, confirms our finding: the binary fraction among FG stars is about 
~12% to be compared with ~1% of SG stars. This is an evidence that SG stars formed in a denser 
environment, where infant mortality of binary systems was particularly efficient. 

Subject headings: binaries: general — globular clusters: general — stars: abundances — stars: Pop- 
ulation II 



1. INTRODUCTION 

About 25 years ago, Renzini & Buzzoni (1986) indi- 
cated globular clusters (GCs) as the best examples in na- 
ture of a simple stellar population, being defined as an as- 
sembly of coeval, initially chemically homogeneous, sin- 
gle stars. However, recent photometric and spectroscopic 
studies significantly challenged this traditional paradigm, 
changing our perspective. In a recent survey, focusing 
on 19 galactic GCs, Carretta et al. (2009b, hereafter 
C09) have convincingly shown that each GC is composed 
of multiple stellar generations, with the coexistence of 
both first (FG, with low Na/Al/N and high O/Mg/C) 
and second generation (SG) stars (high Na/Al/N and 
low O/Mg/C). Briefly, according to Na abundance and 
[O/Na] ratios, GC members were classified as FG (Pri- 
mordial (P) stars) and SG stars, the last one being fur- 
ther divided into Intermediate (I) and Extreme (E; see 
Carretta et al. 2009b for a detailed description on PIE 
classification). The main results found by Carretta and 
coworkers can be summarized as follows: (1) the P popu- 
lation is present in all GCs (at the ^30% level), (2) the I 
population constitutes the bulk (50%- 70%) of GC stars, 
and (3) the E component is instead not present in all 
GCs. 

The chemical composition of SG stars indicates that 
they are formed from FG stars within a limited mass 
range; hence, in order to reproduce the currently ob- 
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served ratio between FG and SG stars (about 1/3 FG 
versus 2/3 SG), the original cluster population should 
have been much larger than the current value. In a com- 
panion paper, Carretta et al. (2010) indeed concluded 
that during the early epochs of dynamical evolution, a 
proto-GC should have lost up to ^90% (or even more, 
see, e.g., Gratton & Carretta 2010) of its P population. 
As a consequence they suggested that the P stars in GCs 
might be the major building blocks of the stellar halo. 

In general, the basic idea expressed in Carretta et al. 
(2010) is that a GC, a few 10 7 yr old, should appear 
as a rather compact aggregate of stars immersed in a 
loose, wide association and in an even larger expanding 
gas cloud (note that such objects are indeed observed in 
extragalactic systems, and especially in the interacting 
ones, see, e.g., Vinko et al. 2009). 

For the subsequent evolution, following the approach 
by D'Ercole et al. (2008), Carretta et al. suggested that 
SG stars formed in the central cluster regions, where 
ejecta of the FG component are channeled into a cool- 
ing flow. A small fraction of P stars are trapped in the 
dense, central cluster regions formed by SG stars: this is 
the typical GC that we can observe nowadays. 

This proposed scenario is still qualitative, given the 
lack of more stringent observational evidences to quan- 
titatively constrain such indications. Along with uncer- 
tainties on cluster formation and early evolution, the na- 
ture of polluter stars responsible for the peculiar chemi- 
cal pattern observed in GCs is also still unclear. In fact, 
the presence of Na— O (and also Mg— Al) anticorrelations 
also among the main-sequence TO stars (Gratton et al. 
2001, Cohen et al. 2002, Carretta et al. 2004, Pasquini 
et al. 2005, Lind et al. 2009, D'Orazi et al. 2010) defini- 
tively demonstrated that a previous generation of stars, 
in which the complete CNO cycle has operated, orig- 
inates chemical (anti) correlations found in GCs. The 
nature of the polluters is still debated, the main candi- 
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dates being intermediate-mass asymptotic giant branch 
(IM— AGB; 5— 8Mq) stars in hot bottom burning phase 
(Ventura & D'Antona 2009), or Fast Rotating Massive 
stars (FRMs; e.g., Decressin et al. 2007), or massive 
binaries (de Mink et al. 2009). 

A complementary approach to proton-capture element 
abundance determination, is the study of the s-process 
elements, like e.g., barium (see, e.g., Smith 2008). This 
information can allow to address both (1) the formation 
and early evolution mechanisms in GCs and (2) the origin 
of the massive, previous generation stars which produced 
the ejecta for the intracluster enrichment. 

The derivation of Ba abundances for a very large sam- 
ple of GC stars allows discovering the presence of the 
so-called Ba stars. Spectra of Ba stars are character- 
ized by unusually strong CH, CN, and heavy elements 
(e.g., Ban, Srli) features. It is now well assessed that 
Ba stars 6 are long-period, single-lined spectroscopic bi- 
nary systems (McClure 1989): the nucleosynthetic pat- 
tern observed can be attributed to now unseen compan- 
ions, which at earlier epoch were low-mass (<3M Q ) AGB 
stars that transferred processed material to the surviv- 
ing visible stars (Luck & Bond 1991). The binary system 
should originally have had a wide enough separation to 
allow evolution of the primary up to the AGB (for a dis- 
cussion on Ba stars in GCs, see Gratton et al. 2004). 

The presence and the fraction of Ba stars within GCs 
indirectly provide hints on binarity properties in dense 
stellar systems, shedding light on formation and evolu- 
tion mechanisms. In this sense, given the quite limited 
number of known Ba stars in stellar clusters, large sam- 
ples are crucial to derive statistically significant infor- 
mation. Moreover, the simultaneous determinations for 
the very same stars of Na,0 and Ba abundances offer the 
possibility to unveil possible differences (or analogies) be- 
tween binary fraction in FG and SG stars. Interestingly, 
determinations of the binary fractions in different stellar 
generations allow to infer key information on density con- 
ditions of their formation environment. In fact, it is now 
well known that the binary incidence is directly related 
to the environmental density (e.g., Lada & Lada 2003): 
this is a basic constraint for hydrodynamical simulations 
of GC formation and early evolution. 

As far as point (2) is concerned, our work aimed at 
determining Ba variation inside each GC as signature of 
low-mass AGB stars contribution to pollution scenario. 
In fact, for stars with masses <3 M Q , the third dredge-up 
becomes very efficient resulting in an s-process element 
enhancement. 

In this Letter we present our results on Ba determina- 
tion for a sample of 15 GCs, for a grand total of more 
than 1200 stars. The huge database allowed us to obtain 
very interesting results minimizing the impact of the ob- 
servational uncertainties. This is by far the largest survey 
of Ba abundances in GC stars available to date. 

2. SAMPLE AND ANALYSIS 

We analyzed high-quality, high-resolution 

FLAMES/Giraffe spectra (Pasquini et al. 2002) 
for a sample of 1205 red giants in 15 Galactic GCs, 

6 Ba stars can be red giants, but also main sequence or sub-giant 
stars; the fact that to date all the Ba stars observed in GCs are 
giants simply reflects a selection effect. 
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Fig. 1. — Run of average [Ba/Fe] with [Fe/H] for our GCs. Empty 
and filled squares are for disk/bulge GC and inner halo GCs; the 
filled triangle is for the only outer halo cluster NGC 1904. Error 
bars for [Ba/Fe] and [Fe/H] are similar to the symbol size. Gray 
symbols are field stars from the literature (Fulbright 2000, Burris 
et al. 2000, Mashonkina k, Gehren 2001, Mashonkina et al. 2003). 
Some GCs are labeled (see the text). 



in order to derive Ba abundances. Details on target 
selections, observations, data reduction, along with 
p-capture elements abundances are given in C09. 

Using the Kurucz (1993) model atmospheres with the 
ROSA abundance code (Gratton 1988), we obtained Ba 
abundances by measuring equivalent widths (EWs) for 
the only Ba II feature covered by our spectral setup 
(HR13), namely the line at 6141 A. 

Ba lines present hyperfine structure occurring in odd 
isotopes ( 135 Ba and 137 Ba); moreover the lines due to the 
five Ba isotopes have small isotopic wavelength shifts. 
However, both these effects are significant only for Ba II 
resonance lines, i.e., 4554 and 4934 A, while the 6141 
A line is scarcely affected (see, e.g., Sneden et al. 1996). 
Moreover, the blend with the Fe I line at 6141.73 A has 
a negligible impact at the low metallicity of GCs and can 
be ignored. 

Iron abundances were presented in Carretta et al. 
(2009a), while stellar parameters (T c g, log g, £) were 
previously derived by C09, to which we refer also for 
details on error estimates. Here, we just mention that, 
since we aimed at discovering possible star-to-star vari- 
ations within each GC, systematic errors have a rather 
small effect. Finally, since this Ba line is near the satura- 
tion regime (with rather high EW values) , the dominant 
source of (internal) uncertainty is given by the adopted 
microturbulences, the final values of effective tempera- 
tures and gravity having only a small impact. 

3. RESULTS AND DISCUSSION 

In Table [U we summarize our results for Ba abun- 
dances in our sample GCs; along with metallicity [Fe/H] 
as derived by Carretta et al. (2009a), we list the abso- 
lute magnitude of the clusters (from Harris 1996), and 
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Fig. 2. — [Na/O] distributions for Ba-rich (red) and Ba-poor 
(black) stars. (A color version of this figure is available in the 
online journal.) 



TABLE 1 
Ba Abundances for Our Program Clusters. 



Cluster 



Mv 



[Fe/H] 
±err 



[Ba/Fc] 
±err 



NGC 104 


-9.42 


-0.743±0.010 


0.150±0.014 


0.150 


110 


NGC 288 


-6.74 


-1.219±0.013 


0.450±0.021 


0.183 


75 


NGC 1904 


-7.86 


-1.544±0.011 


0.236±0.026 


0.185 


49 


NGC 3201 


-7.46 


-1.495±0.014 


0.236±0.026 


0.267 


104 


NGC 4590 


-7.35 


-2.227±0.020 


-0.128±0.033 


0.316 


93 


NGC 5904 


-8.81 


-1.346±0.006 


0.095±0.021 


0.210 


108 


NGC 6121 


-7.20 


-1.200±0.007 


0.502±0.013 


0.118 


80 


NGC 6171 


-7.13 


-1.065±0.020 


0.612±0.041 


0.215 


27 


NGC 6254 


-7.48 


-1.556±0.014 


0.156±0.030 


0.305 


105 


NGC 6397 


-6.63 


-1.993±0.011 


-0.024±0.030 


0.306 


103 


NGC 6752 


-7.73 


-1.561±0.011 


0.277±0.017 


0.176 


106 


NGC 6809 


-7.55 


-1.967±0.012 


0.166±0.027 


0.284 


112 


NGC 6838 


-5.60 


-0.808±0.010 


0.269±0.027 


0.143 


29 


NGC 7078 


-9.17 


-2.341±0.017 


0.189±0.055 


0.412 


56 


NGC 7099 


-7.43 


-2.359±0.015 


-0.218±0.021 


0.142 


18 



the mean Ba abundance for each GC with the corre- 
sponding error and the number of stars. In general, the 
rms for individual stars from the average [Ba/Fe] ratios, 
0.15-0.30 dex, are of the same order of magnitude of the 
expected errors (about 0.25-0.30 dex), pointing out that 
there is no star-to-star variation in Ba content, save for a 
few exceptional stars (see discussion in the following sec- 
tion). The only exception is NGC 7078, which shows the 
largest standard deviation from the mean, i.e., 0.41 dex. 
For this cluster, we derived [Ba/Fe]=0.189±0.055 to be 
compared with [Ba/Fe]=0.12±0.04 (rms=0.21) inferred 
by Sneden et al. (2000) from a sample of 31 giants. The 
quite high [Ba/Fe] ratio along with a significant disper- 
sion for NGC 7078 confirm the previous study by Sneden 
et al. (2000; see also Armosky et al. 1994 for the first 
suggestion of this trend). However, the [Ba/Eu] ratio 
derived by Sneden et al. (1997, 2000) indicates that the 



heavy element abundances for this GC are compatible 
with a pure r-process pattern: the bimodal distribution 
for Ba and Eu reveals an r-process origin both for Ba 
and Eu. For the other GCs for which Ba abundances are 
available in literature, our estimates are in good agree- 
ment with the previous ones. Among the others, for 
the intermediate-metallicity GC NGC 6121 we obtained 
[Ba/Fe]=0.502±0.013, which is very close to the value 
derived by Ivans et al. (1999; [Ba/Fe] =0.6 dex). Also, 
we confirm that NGC 6121 shows an overabundance in s- 
process elements with respect to its twin NGC 5904 (for 
which we derived [Ba/Fe] =0.095±0. 021), providing sup- 
port to previous estimates (see, e.g., Ivans et al. 1999, 
2001). These comparisons with previous determinations 
prove that no major systematic uncertainties affect our 
analysis. 

In Figure [TJ we show the run of [Ba/Fe] as function 
of metallicity for our 15 GCs, with different symbols 
separating disk/bulge, inner halo, and outer halo GCs 7 . 
There is not a clear separation between the different 
Galactic components, with a large scatter in [Ba/Fc] ra- 
tio at any given metallicity. On the other hand, the well- 
known raising trend of s-process with [Fe/H] is easily 
recognizable. This can be attributed to the AGB stars 
contribution to the Ba production, which becomes dom- 
inant at metallicity near [Fe/H]« —1 (see Busso et al. 
1999). 

None of the clusters present a correlation between Na 
and Ba abundances (see also James et al. 2004), indi- 
cating that there is no significant contribution from low- 
mass AGB stars. We could reach the same conclusion by 
"considering the cumulative distributions for Ba-rich and 
Ba-poor stars: 8 we found that the two [Na/O] distribu- 
tions are almost the same, with a Kolmogorov— Smirnov 
test returning a probability of 99.8% that they are ex- 
tracted from the same population (see Figure [2]). Also 
by dividing our sample stars in P, I, E (C09), the three 
cumulative curves for [Ba/Fe] ratios are very similar: 
within ^1.5<t the [Ba/Fe] values for FG and SG stars 
are the same. 

3.1. Ba stars 

Despite the quite large uncertainties («0.25 dex), 
mainly due to the fact that our spectra were not acquired 
with the specific purpose of performing heavy element 
_abundance analysis, the large statistics lead to the dis- 
covery of five Ba stars. We defined a Ba star if [Ba/Fe] 
is more than 3.5er above the average of the other stars. 
As an example, we compare in Figure [3] the spectra of 
two stars (# 28903 and # 28881, both in NGC 6254) 
which have the same T e g, log g and [Fe/H]; the first is 
Ba rich, while the latter has normal Ba value. We found 
that the EWs of Ba II lines differ by more than a factor 
of two. Given the similarities in atmospheric parameters, 
different EW strengths must correspond to different Ba 
content. For our five Ba stars , other s-process elements 
were determined, whenever possible, as reported in Ta- 
ble [2j Discarding elements characterized by features in- 

7 GCs are divided according to a combination of positional and 
kinematics criteria. A detailed description of the followed proce- 
dure is given in Carretta et al. (2010). 

8 For each GC, Ba-rich and Ba-poor stars are defined as stars 
with [Ba/Fe] larger and smaller, respectively, than the median 
value. 
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trinsically weak in our stars, we derived abundances for 
La II and Nd II. Taking only lines with EWs larger than 
10 ml, we used the 5769.06 A, 5797.57 A, and 6390.48 
A lines for La II, while the 5740.86 A and 5804.00 A ones 
for Nd II. 

Over the whole sample of 1205 stars, we derived a frac- 
tion of Ba stars of ^0.4%, confirming previous sugges- 
tions that the fraction of Ba stars in dense stellar systems 
like GCs is significantly lower than in the field (~2%, 
Luck & Bond 1991). Furthermore, Gratton et al. (2004) 
argued that all known CH and Ba stars in GCs are in 
low concentration clusters: this is probably because in 
highly concentrated clusters, the possible progenitor sys- 
tems are destroyed on a timescale shorter than the evo- 
lution of intermediate low-mass stars (Giersz & Spurzem 
2000, see also Gratton et al. 2004). However, we found 
no relationship between Ba stars and cluster concentra- 
tion, with Ba stars present in both NGC 288 and 47Tuc 
(NGC 104), which have central concentration parameters 
of c=0.96 and c=2.03, respectively (see Harris 1996). In 
addition, we did not find any correlation with metallicity 
and absolute magnitude. 

In Table [2j along with heavy elements, we report for 
our Ba stars also the abundances for Na and O mea- 
sured by C09. Four out of the five Ba stars detected in 
our sample belong to the P population. This means that 
if we compute the fraction of Ba stars only considering 
the P stars, the percentage is ~ 2%, the same number 
observed in field stars. This interesting finding supports 
our scenario for GC formation, where the SG stars were 
born in a denser environment (generated by the cool- 
ing flow), while the primordial population formed as an 
extended and loose association. The environment condi- 
tions could have influenced the wide binaries population, 
since in high-density conditions the high rate of colli- 
sions significantly contribute to the destruction of such 
systems, while in sparse environment binaries can survive 
in a larger number. 

As an independent check, to confirm the suggestion 
from Ba stars, we also analyzed radial velocity variations 
for a hundred giant stars in NGC 6121, looking for spec- 
troscopic binaries (more details will be given in a sepa- 
rate paper, S. Lucatello et al., 2010, in preparation). We 
chose stars for which at least two radial velocity deter- 
minations were available, either within our observations 
program (time span of At=3 months) or combining our 
data and the observations carried out by Marino et al. 
(2008; At=2 years). We discovered five binary stars on a 
grand total of 102, which means a fraction of ~5%: our 
result is in good agreement with the binary fraction for 
this GC found by A. P. Milone et al. (2010, in prepara- 
tion), i.e., 6.0%±0.4% (by considering only systems with 
a mass ratio g=Ml/M2 >0.5, see Milone et al. 2008 for 
further details on this point.) 

Considering the O and Na content for these binaries, 
we found that four out of five belong to the P component 
(with low Na content). We performed a statistical test 
and the probability that this finding can be obtained by 
chance is only of ~5%. Assuming a fraction of 34%±6% 
for FG and 66%±8% for SG stars 9 , we hence derived 

9 The 102 stars of our sample were divided according to the 
Na abundance, with FG having [Na/Fe]<0.28 dex (see the Intro- 
duction). Our estimate is in good agreement with the values of 
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Fig. 3. — Comparison for the spectra of Ba star 28903 (in 
NGC 6254) with a normal star with very similar atmospheric pa- 
rameters (star 28881). Note the strong differences in Ba II line 
strengths. 



a value of ~12% for the frequency of binaries among P 
stars, and only ~1% for SG stars. 

Finally, we investigated whether there is a trend be- 
tween the fraction of FG stars within a GC and the 
number of Blue Stragglers (BSS), since BSS seem to have 
predominantly an origin as P binaries, with only a small 
percentage formed through different mechanisms, e.g., 
collisions (see, e.g., Moretti et al. 2008). By comparing 
the fraction of P stars in 14 GCs in common between our 
previous paper (C09) and the one by Moretti et al., we 
indeed find that there is a positive correlation between 
the number of FG stars and the BSS fraction (the linear 
correlation coefficient has a significance level of 97.5%). 
This is a further confirmation of what derived from both 
Ba stars and spectroscopic binaries in NGC 6121. 

Our result is yet another observational evidence, along 
with metallicity distribution, mass and location, of the 
similarities between the primordial component of GCs 
and the field stars. In a forthcoming paper (R. G. Grat- 
ton et al., 2010, in preparation), we will widely discuss 
this issue. Since second generation stars are nowadays 
dominant in GCs, our result could explain why binary 
fractions in GCs are smaller compared to the field (e.g., 
Pryor et al. 1989; Hut et al. 1992; Cote et al. 1996). 
Note that at variance with the field, different processes 
determine the relative frequency of binary systems in 
stellar clusters; in GCs binaries are continuously formed 
and destroyed during the cluster evolution as the result of 
collisional interactions between binaries and single stars 
(see, e.g., Sollima et al. 2009). However, the cluster re- 
gions sampled by FLAMES are generally quite far from 
the center, where dynamical interactions are much more 
frequent. This might explain why we can observe these 

30%±5% and 70%±8% found by C09, respectively for FG and SG 
stars. The main difference between our computation and the one 
by C09 is that we divided FG and SG using only Na abundances, 
while C09 considered both Na and O contents. 
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differences between stars belonging to different stellar 
generations some 10-13 Gyr after GC formation. 

4. SUMMARY AND FUTURE PERSPECTIVES 

In this Letter, we report our results on Ba abundance 
determination for a sample of more than 1200 giants in 
15 Galactic GCs. We provide the largest, homogeneous 
analysis of this kind available in literature to date. 

We found that there is no star-to-star variation in Ba 
content within each cluster, with the standard devia- 
tion from the mean compatible with the expected un- 
certainties (the only exception being the low-metallicity 
NGC 7078, for which, however, a pure r-process pat- 
tern has been proposed by several authors). Our results 
suggest that there is no detectable contribution to intra- 
cluster pollution from low-mass AGB stars, whose prod- 
ucts are s-process elements enhanced due to the efficiency 
of the third dredge-up. 

We discovered the presence of five Ba stars, confirm- 
ing that the fraction of these originally quite wide binary 
systems is lower in GCs than in the field. We also found 
that four out of the five Ba stars belong to the P popula- 
tion: this likely reflects the different environmental con- 
ditions where the two stellar generations formed. This 
new result supports the formation/evolution scenario for 
GC in which the SG stars were born at the cluster cen- 



ter, where the infant mortality of long-period binaries is 
significantly higher. 

To have a further validation of this piece of evidence, 
we also searched for binary stars in NGC 6121: for this 
cluster, we determined rather large radial velocity varia- 
tions for five stars on a total of ~100. Also in this case we 
obtained that only one binary belongs to SG stars, and 
that the fraction of binary stars among the P component 
reaches a value of ~12%, to be compared with the small 
percentage in SG, i.e., ~1%. 

In the future, we will focus on the determination of 
radial velocity variations among a large sample of GCs 
and of stars within each GC. We aim at inferring binarity 
properties in the FG and SG of GCs, as excellent tracer 
of the environmental conditions where GCs formed and 
evolved. 
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D'ORAZI ET AL. 



TABLE 2 
Abundances for Ba-rich stars. 



Star 



cluster 



[O/Fc] 



[Na/Fc] 



[Ban/Fc] 



[Lan/Fc] 



[Ndn/Fc] 



30952 NGC 104 

200095 NGC 288 

28903 NGC 6254 

608024 NGC 6397 

38291 NGC 6752 



0.336±0.046 
-0.045±0.088 
0.428±0.036 
0.387±0.048 
0.253±0.043 



0.343±0.023 
0.079±0.053 
0.203±0.031 

0.126±0.038 



0.508±0.150 
1.115±0.210 
1.464±0.160 
1.040±0.344 
1.261±0.182 



1.209±0.044 
0.281±0.150 



1.035±0.030 



